Compensated magnetism by design in double perovskite oxides 
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Taking into account Goodenough's superexchange rules, including both full structural relaxation 
and spin-orbit coupling, and checking strong correlation effects, we look for compensated half metals 
within the class of oxide double perovskites materials. Identifying likely half metallic (or half 
semimetallic) antiferromagnets, the full complications including orbital magnetism are included 
in order to arrive at realistic predictions of designed magnetic compounds with (near) vanishing 
net moment. After sorting through several candidates that have not been considered previously, 
two materials, K^MnRhOe and I^CrWOe, remain as viable candidates. An important factor is 
obtaining compounds either with very small induced orbital moment (helped by closed subshells) 
or with an orbital moment that compensates the spin-orbit driven degradation of half metallic 
character. While thermodynamic stability of these materials cannot be ensured, the development 
of layer-by-layer oxide deposition techniques does not require that materials be thermodynamically 
stable to be synthesized. 

PACS numbers: 



I. INTRODUCTION 

Double perovskites have drawn considerable at- 
tention in the field of spintronics, particularly 
since the discovery-i of colossal magnetoresistance 
in S^FeMoOe that is believed to be a half metallic 
(HM) ferromagnet (FM). The HM property is highly 
desired in spintronics applications, where spin cur- 
rents are utilized as well as charge currents. The 
double perovskite family of oxides has been one of 
the most popular classes within which to look for 
HMFMs.- 

Recently, with wider application of deposition 
techniques leading to better materials and improved 
understanding, layer by layer deposition of per- 
ovskite oxides has led to unusual and potentially 
useful physical properties. 3 Most such depositions 
of perovskites have used the (001) growth direction, 
resulting in a wide variety of oxide heterostructures 
with tetragonal symmetry. Double perovskite ma- 
terials can be considered as single unit cell, multi- 
layered perovskite structures that have been grown 
along the (111) direction, and such growth is a 
promising direction for synthesizing new members 
of this class. 

Although HMFM character will provide new (and 
highly sought after) magnetoelectronic properties, 
there is a subclass that is more exotic. When half 
metallicity occurs and in addition the moments com- 
pensate, a "compensated half metal" (historically 
called a half metallic antiferromagnet 4 ) arises. This 
state has magnetic order but vanishing macroscopic 
spin magnetization, it could support a new type of 
and one can imagine numerous 



superconductivity, 5 



possible applications for a half metal that is rela- 
tively impervious to external magnetic fields. How- 
ever, this spin-only picture is an idealization; spin- 
orbit coupling (SOC) couples the two spin directions 
and thereby destroys the precise spin compensation, 
and in addition generates spin-induced orbital mo- 
ments. Thus a fully compensated half metal with 
vanishing macroscopic magnetization is an idealiza- 
tion, and the focus must be on finding realizations 
that are as close as possible to full compensation. 

Interest has been rekindled recently in the effects 
of spin-orbit coupling in oxides because, unlike the 
degradation of half metallicity that is commonly dis- 
cussed and expected, it can also lead to more exotic 
effects in the spin- and orbital-magnetism, especially 
in d 1 or d 5 ions in ti g subshells where large or- 
bital moments may arise^ A dramatic possibility is 
the promotion of compensating spin and orbital mo- 
ments. In Ba2NaOs06, for example, calculations 7 
indicate that the d 1 spin moment of Os is compen- 
sated by the ti g L=l orbital moment induced by 
SOC, and it is the partial orbital quenching by the 
environment that destroys the compensation of the 
moment. In Sr 2 CrOs06, neglect of SOC leads to a 
semimetallic compensated half metal, which is de- 
stroyed by the large SOC in this system^ In isova- 
lent Si^CrRuOe, where weaker SOC arises on Ru 
compared with Os, the reduced spin and induced 
orbital moments conspire to nearly cancel, giving a 
real possibility for an effectively compensated half 
metal. In correlated insulators, strong SOC can lead 
to much more complex magnetic coupling, thereby 
complicating the resulting magnetic order— and re- 
ducing the likelihood of achieving a compensated 



half metal. 

Half-metallic antiferromagnetic (HMAF) ma- 
terials have not been obtained experimentally 
yet. Several theoretical efforts have been car- 
ried out in predicting half-metallic antifcrromag- 
netism, but experimentally, this state remains 
elusive. HMAFs have been predicted for vari- 
ous double perovskiteS ) 1 ^ 11 ' 12 ! 13 tetrahedrally co- 
ordinated compounds 1 ^ and for Heusler-structured 
materials i 15 ' 16 i 17 Design of materials with desired 
properties is a long-standing hope that is gaining 
momentum^ and the considerable experience that 
has been accumulated in oxide double perovskitcs 
provides guidelines to focus design efforts. 

In this paper we build on previous studies, and 
extend the earlier work by considering all of factors: 
ion size, structural relaxation, and SOC. Our gen- 
eral strategy is: 

i) we focus on the double perovskite structure, with 
the chemical formula A2MM'Og, where M and M' 
are transition metal ions lying on a rocksalt type 
sublattice. An ordered double perovskite is more 
likely to occur if M and M' differ in ionic radius, 
and differ in formal charge. 

ii) we use ab initio calculations to relax the struc- 
ture, including volume optimization, c/a optimiza- 
tion, and internal coordinate relaxation, thereby ob- 
taining structures that are dynamically stable, and 
energetically metastable if not stable. 

iii) we include SOC, which mixes the spin moment 
with an orbital moment, and also analyze the effects 
of correlation corrections by using the LSDA+U 
method ) 19 ' 20 to check the dependence of the mag- 
netic and electronic structure properties on U (the 
on-site Coulomb repulsion). The objective is to 
determine to what degree the spin-compensation 
is maintained, and to understand the interplay of 
structural changes, SOC, and effects of electronic 
correlation. 

We use the Goodenough-Kanamori-Anderson 
(GKA) rules for superexchange 2 ^ for the 180° metal- 
oxygen-metal bond occurring in perovskites to guide 
our choices. According to these authors, and the ex- 
pected moments, only a few d™-d m electronic con- 
figurations that would lead to a spin-compensated 
situation are likely. In addition, we aim to maximize 
the likelihood of obtaining an ordered perovskite by 
mixing a 3d row element with another metal cation 
from the 4d or 5d row, so that the size mismatch 
will make chemical ordering more likely. According 
to GKA, the possible antiferromagnetically coupled 
electronic configurations are: d 2 -d 2 , d 2 -d 4 (LS), d 3 - 
d 3 , d 8 -d 8 and other combinations that include in all 
cases high spin (HS) cations: d 4 -d 4 , d 5 -d 5 , d 4 -d 6 , 
d 6 -d 6 , d 7 -d 7 . These latter options have not been 
pursued because it is difficult to stabilize HS states 



with those large moments for the 4d and 5d ele- 
ments, which usually take on a higher valence state 
than the 3d element, and hence the relatively strong 
crystal field will favor low-spin states. Particularly 
interesting is the d 3 -d 3 configuration when SOC is 
taken into account. Since it deals with a completely 
filled t% g band, the effective angular momentum will 
be very small, although there is also the possibility 
of insulating states. We will discuss this issue later 
in the paper. 



II. COMPUTATIONAL DETAILS 

Electronic structure calculations were performed 
within density functional theory^ 2 - using wien2k^ 
which utilizes an augmented plane wave plus local 
orbitals (APW+lo) 24 method to solve the Kohn- 
Sham equations. This method uses an all-electron, 
full-potential scheme that makes no shape approx- 
imation to the potential or the electron density. 
The exchange-correlation potential used was the 
Perdew and Wang version of the local density 
approximation^ and strong correlation effects were 
introduced by means of the LSDA+U schem o 26 ' 27 
including an on-site U (Coulomb repulsion) and J 
(Hund's rule exchange) for the metal d states. We 
have used option nldau= 1 in wien2k, i.e. the so- 
called "fully-localized limit" . Structural optimiza- 
tions (unless otherwise stated) were performed using 
the Perdew-Burke-Ernzcrhof version of the general- 
ized gradient approximation. 28 Spin-orbit coupling 
was introduced using the scalar relativistic approxi- 
mation. All the calculations were converged with re- 
spect to all the parameters used, to the precision nec- 
essary to support our calculations (converged forces 
and total energy differences), up to H mt K max = 7 
and an 8 x 8 x 8 k-mesh. 



III. REVISITING PREVIOUSLY FOUND 
MATERIALS 

A. La 2 VCu0 6 : cP-d 9 

It had been reported previously that this material 
is a HMAF ) 10 ' 29 but previous studies did not con- 
sider the structural relaxation that we have now car- 
ried out, and that is expected to be important in this 
strongly Jahn- Teller active system. At the LSDA+U 
level, this compound converges to a HMAF, which 
remains when one relaxes the structure. We have 
done a relaxation within a tetragonal symmetry (but 
the ratio c/a ~ 1.00, to 1% accuracy, as expected), 
obtaining a= 3.64 A. The internal structure relaxes 
to accommodate the Jahn- Teller distortion of the 
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Cu 2+ and V 4+ cations. Depending on the value of 
U chosen (the same U was chosen for both cations in 
our calculations), this distortion varies, being larger 
for a larger U. For relatively small values of U (ef- 
fective U of only 4 eV) there is a large Jahn- Teller 
distortion, the oxygen octahedra around the Cu 2+ 
cation become elongated by about 4% to accommo- 
date the d 9 cation and around V 4+ , the octahedron 
is shortened along one axis by about 4% to acco- 
modate a d 1 cation. The system becomes a Mott 
insulator for U= 4.8 eV and J= 0.7 eV. The elec- 
tronic structure is presented in Fig. [T] for U= 3.5 
eV, where the system is still half-metallic. 

Figure [T] shows that half-metallicity occurs with 
the Fermi level falling within a 2 eV wide band of 
Cu e g and V t2 S character, even when degeneracy is 
broken by the Jahn- Teller distortion. The bottom 
of the band (a two band complex) is strongly Cu d 
in character, while the V d character dominates the 
upper part of the bands. 

We expect a small influence of SOC because only 
3d elements are involved. The full moment com- 
pensation that occurs within LSDA+U will be kept 
when including SOC. This d 9 -d 1 configuration has 
small magnetic coupling according to GKA rules 
(AF according to our calculations), and is spin- 
compensated according to our results. The orbital 
angular momentum of the Cu atoms is 0.1 /is, 
but is negligible for V and the electronic structure 
is not affected by the introduction of SOC. Spin- 
compensation remains unaltered when SOC is in- 
troduced, thus the total moment of the system is 
approximately 0.1 //s/f.u., quite close to the desired 
result. 



B. Sr 2 NiOs0 6 : d 8 -d 2 

This material has been synthesized 30 recently and 
also has been studied using ab initio methods. 31 
Experimental results were interpreted in terms of 
a competition between FM and AF interactions, 
caused by the d 8 -d 2 (Ni 2+ -Os 6+ ) electronic struc- 
ture. The closed subshell (d 8 ), open subshell (d 2 ) 
should favor half metallicity. According to GKA 
rules, the e 2 -0-e°, FM superexchange should be the 
strongest in this system. However, the magnetic sus- 
ceptibility data shows both a positive Qcw (show- 
ing that FM is dominant at high temperature) and 
an AF downturn below 50 K, indicating that AF 
correlations exist at low temperature. Our ab ini- 
tio calculations show the FM and AF configurations 
are similar in total energy, but FM is distinctly the 
ground state by 20 meV/metal atom (using the ex- 
perimental structure). Also, the spin-compensation 
that results within LSDA+U is degraded seriously 
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FIG. 1: Density of states plots of the tetragonal dou- 
ble perovskite compound La2VCu06 in a HMAF state 
within an LSDA+U scheme with U= 3.5 eV for both 
V and Cu. Upper (lower) panel shows up (down) spin 
channel. The spin moment is precisely compensated. 

by SOC, with the large orbital angular momen- 
tum leading to a total moment of approximately 1 
/xs/f.u. This value arises from orbital moments of 
0.5 /jb for Os and 0.2 jig for Ni, plus the 0.2 /zg 
from the total spin moment of the system arising 
from SOC-induced spin mixing. Hence, the condi- 
tion of a small magnetic moment is not obtained for 
this compound. The d 2 -d 8 electronic structure has 
the drawback that, when an antialigned spin con- 
figuration is set, the orbital angular momenta align 
and add up rather than (partially) canceling (Hund's 
third rule for the atoms). 

We have performed a structural analysis similar 
to the one described for La2VCu06 with the same 
preconditions. Comparing with experimental struc- 
tural data^ our calculated volume is quite good 
(1% larger than experiment) and the c/a ratio is 1% 
smaller than the experimental value. 

TABLE I: Calculated versus experimental structural pa- 
rameters and selected bond lengths for Sr2NiOs06 show- 
ing agreement to within 1% accuracy in the lattice con- 
stants. 
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Experimental 30 


This work 


a 




5.49 


5.53 


c 




7.99 


7.94 


Ni-O distance (x 


4) 


2.02 
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Ni-O distance (x 


2) 


2.04 


2.03 


Os-O distance (x 


4) 


1.91 


1.92 


Os-O distance (x 
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1.96 


1.94 
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FIG. 2: Density of states plots of the tetragonal dou- 
ble perovskite Sr2NiOs06 in a HMAF state within an 
LSDA+U scheme with U= 5 eV for Ni and 3 eV for 
Os, including spin-orbit coupling. Upper (lower) panel 
shows up (down) spin channel. 

For large values of U this compound becomes a 
Mott insulator. Half-metallicity occurs only for U 
smaller than 3.5 eV for Os (with J= 0.7 eV) and 
values of U not larger than 7 eV for Ni. Realistic val- 
ues of U and J would be close to the metal-insulator 
transition. No experimental information on its con- 
duction properties is available so far. For smaller 
values of U, we obtain conduction coming from the 
Os majority spin channel (the unfilled subshell), as 
can be seen in Fig. [2] 

IV. NEW FINDINGS 

Two materials with the desired properties have 
been clearly identified. Both possess a d 3 -d 3 elec- 
tronic configuration so a high-spin, magnetic band 
insulator situation is a possibility. This particu- 
lar electronic configuration has several advantages, 
that we will discuss, over other possible electronic 
configurations that will promote a spin-compensated 
state, and they have AF coupling by the GKA rules, 
so they become plausible candidates for being spin- 
compensated half metals or half semimctals. In ad- 
dition, we have analyzed a compound with a d 2 -d 2 
configuration, and we discuss distinctions. 

A. Sr 2 VRe0 6 : d 2 -d 2 

This compound is calculated to have a d 2 -d 2 elec- 
tronic configuration. The unoccupied t^g shells lead 



to a distortion from cubic symmetry. Relaxing the 
structure within tetragonal symmetry, the resulting 
c/a ratio is 1.02, together with a sizable elongation 
of the oxygen cages to accommodate the d 2 elec- 
tronic state (3% difference in the V-0 distance and 
a smaller 1% difference for the Re-0 distance). A 
summary of the structural parameters obtained is 
given in Table [Til 

TABLE II: Calculated structural parameters for 
Sr2VRe06 in the space group 14/mmm (no. 139). a— 
5.52 A( v / 2 x 3.90 A), c=7.96 A(2 x 3.98 A). 



Atom 


Wyckoff position 


Atomic parameter 


Sr 


4d 


(0.5,0,0.25) 


V 


2a 


(0,0,0) 


Re 


2b 


(0,0,0.5) 


Ol 


8h 


(0.2521,0.2521,0.5) 


02 


4e 


(0,0,0.2494) 



We have studied the electronic structure of this 
yet-unreported compound with both LSDA+U and 
SOC. Due to slight band overlap, there is substantial 
but not total spin-moment compensation, leading to 
a total spin moment of less than 0.1 fiB- However, 
because of the unfilled t2 g shells in both cations, 
large orbital angular moments develop, particularly 
large for Re (0.4 ^b)- The result is a total mag- 
netic moment of approximately 0.5 fiB for U(V)= 7 
eV and U(Re)= 3 eV (J= 0.7 eV for both cations). 
While there is no true gap in the V majority spin 
channel, there is a very large fractional spin polar- 
ization in the density of states at the Fermi level, as 
can be seen in Fig. [3] 

B. K 2 MnRh0 6 : d 3 -d 3 

As we have just seen in the case of S^VReOe, 
having unfilled shells encourages orbital angular mo- 
ments. In that respect, the following two examples 
are materials which have a d 3 -d 3 electronic config- 
uration, so the filled t2 S shells will give very small 
orbital moments to obscure the spin compensation. 

We calculate that this K^MnRhOg compound has 
the configuration Mn 4+ : d 3 - Rh 6+ : d 3 , which cou- 
ples AF. The Rh moment is strongly delocalized on 
the oxygen neighbors (up to 0.3 /ib per oxygen and 
only 1 fis in the Rh muffin-tin sphere), but all in 
all, this compound is near HMAF with Fermi level 
lying very near band edges. The lattice parameter 
a was optimized and atomic positions were relaxed 
within the Fm3m space group (no. 225), and a— 
3.94 A was obtained. The d 3 -d 3 configuration of 
ions with cubic symmetry leaves a cubic structure, 
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C. La 2 CrW0 6 : d 3 -d 3 
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FIG. 3: Density of states plots of the tetragonal dou- 
ble perovskite Sr2VRe06 in a HMAF state within an 
LSDA+U scheme with U= 5 eV for V and 3 eV for Re, 
including spin-orbit coupling. Upper (lower) panel shows 
up (down) spin channel. 



leaving only a small relaxation of the oxygen cages 
around the cations (Rh-0 distance of 2.00 A and 
Mn-0 distance of 1.94 A). 

As mentioned, very small values of the orbital 
angular momentum arise when SOC is introduced: 
0.05 fiB on Rh, 0.01 \ib on Mn. Adding spin and or- 
bital angular momenta, the net moment is only 0.02 
/ib, and it remains very close to compensated half 
metal for values of U of 6 eV for Mn and 4 eV for Rh. 
Even though Fig. 2] looks similar to a zero-gap sit- 
uation, on the verge of a metal-insulator transition, 
bands cross the Fermi level even for larger values of 
U, and (half) semimetal character is maintained. 
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This is another example of a d 3 -d 3 AF coupling 
situation, and again, it is a not-yet-synthesized ma- 
terial, completely ab initio designed. Structural re- 
laxation was carried out allowing tetragonal distor- 
tion. No tetragonal distortion was obtained, and a 
= 4.00 A. The difference in ionic radii leads to a 
Cr-0 distance of 1.97 A and W-0 distance of 2.03 

A. 

Evolution with U and SOC indicates that the lat- 
ter does not seriously degrade the moment compen- 
sation. For large values of U, where the material 
is still a metal, the total moment stays below 0.2 
/zs/f.u. W has an orbital moment of 0.15 fiB and a 
much smaller one, and opposite in an AF configura- 
tion, arises from Cr. Even for large values of U (7 eV 
for Cr and 4 eV for W), the system remains metal- 
lic (half- metallic) and the total magnetic moment, 
including spin and orbital components, does not ex- 
ceed 0.10 /zs/f.u. Metallicity comes from a broad W 
d band that crosses the Fermi level, as can be seen 
in Fig. [SJ An electron pocket is formed in an itin- 
erant W d band and a corresponding La-character 
hole pocket (not shown) appears. 
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FIG. 5: Density of states plots of the tetragonal dou- 
ble perovskite L^CrWOe in a HMAF state within an 
LDA+U scheme with U= 5 eV for Cr and 3.5 eV for 
W, including spin-orbit coupling. Upper (lower) panel 
shows up (down) spin channel. 



FIG. 4: Density of states plots of the tetragonal dou- 
ble perovskite K^MnRhOe in a HMAF state within an 
LSDA+U scheme with U= 5 eV for Mn and 3 eV for 
Rh, including spin-orbit coupling. Upper (lower) panel 
shows up (down) spin channel. 



V. UNPROMISING COMPOUNDS 

On our way to detecting these spin-compensated 
half-metals, we studied other seemingly promising 
double perovskites. We will discuss these only 
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briefly, since they do not show the desired behav- 
ior. 



A. d^d 9 

Even though superexchange rules predict this cou- 
pling to be negligible, we have previously studied a 
system that could be promising. Hence, as an anal- 
ogy with the Jahn- Teller compound I^VCuOg, we 
tried changing each of the 3d elements by a corre- 
sponding 4d element, without changing the electron 
count, i.e. La 2 CuNbOg and La 2 VAg06 which are 
less likely to be Mott insulators. In these cases, Nb 
and Ag respectively are non-magnetic, and the spin- 
compensation is lost. 

B. d 2 -d 4 

This electronic configuration could lead to a spin- 
compensated state, if the d cation is in a low-spin 
state. We have tried Ba2FeCrOg (also with Sr and 
Ca with similar results), which could assume a d 2 -d 4 
electronic configuration. It is indeed the case, but Fe 
is in a HS state (not allowing for spin compensation 
to happen). A similar problem appears with the 
compound K 2 FeRuOg, where also Fe 4+ is also in 
a HS state. This electronic configuration has been 
studied in the pas t 11 ! 12 for several perovskites with 
general formula AA'BB'Og- 

C. d 2 -d 8 

We have found several cases where always the FM 
exchange mediated by the e g electrons is more in- 
tense than the AF coupling between the t2 g elec- 
trons. Starting from Sr 2 NiOs06 compound dis- 
cussed above and substituting Os by Ru to form 
Sr2NiRuC>6, we find that FM coupling remains. This 
is the standard prediction of GKA rules. Apart from 
these, we have also identified that FM coupling is 
more favorable in La2CrNiOg (and also with W and 
Mo on the Cr site). This material is calculated to 
be an interesting FM insulator (of which there are 
relatively few examples), where instead of promot- 
ing the possible Jahn- Teller distorted Ni 3+ cation, 
it leads to the formation of a Cr 2+ -Ni 4+ configura- 
tion. This charge ordered state is thought to occur 
to prevent the system from developing a Jahn- Teller 
distortion when it is close to the itinerant electron 
limit. 32 

Disordered LaCro.5Nio.5O3 is known experimen- 
tally to be insulating, 33 but close to a metal- 
insulator transition (the Ni end of the series LaNiO.3 



is a metal). Our calculations show that the elec- 
tronic structure of the charge ordered material ex- 
plains the insulating state found experimentally, to- 
gether with the prediction of FM coupling between 
the d 2 and d 8 cations, making the system poten- 
tially interesting as another FM insulator. For these 
calculations, we have taken the lattice parameters 
from the experimental with the disordered system 
and relaxed only the internal coordinates. The dif- 
ferent ionic radii lead to different metal-oxygen bond 
lengths (1.90 A for Cr and 1.96 A for Ni). 



D. Antiferrimagnetic insulators 

We have found the following magnetically- 
compensated insulators with double perovskite 
structure: i) La2CrMoOg is isoelectronic to the 
spin-compensated metal La2CrWOg, but introduc- 
ing Mo instead of W opens up a gap at the Fermi 
level. It is the greater itineracy of the 5d element 
that leads to the metallic state that we analyzed 
above, ii) In analogy with Sr 2 CrOs06^ we have 
tried with a bigger alkaline-earth cation. The com- 
pound Ba2CrOs06 is found to be an antifcrromag- 
netic insulator, with configuration d 3 -d 3 . As an iso- 
valent compound, S^CrRuOg is calculated to be an 
antiferromagnetic insulator. 



VI. SUMMARY 

In this paper we have suggested two very nearly 
spin-compensated half metals (or half semimetals): 
K 2 MnRh06 and La2CrW06, both with a formally 
d 3 -d 3 configuration. It is found that, within a dou- 
ble perovskite structure, a d 3 -d 3 electronic configu- 
ration promotes the chances to obtain a nearly spin- 
compensated half- metal, due to several factors, i) 
No competition between FM and AF interactions 
occurs as happens when there is partial occupation 
of e g states in one of the metal cations, ii) Having 
the t 2g shell completely filled leads to a very small 
orbital angular moment, which works even further 
in favor of a very small magnetic moment solution, 
iii) Having less than half-filled d shells in both metal 
ions helps to avoid imbalance in the total moment 
that occurs when spin-orbit coupling is important, 
specially for 4d and 5d elements. Having both cation 
d shells less than half filled means that the orbital 
angular momenta of both cations will also orient an- 
tiparallel to the spin, preventing the system from 
having a larger net magnetic moment, iv) The large 
S=3/2 spins make our calculations more realistic; 
ordering of low spin ions is often more affected by 
quantum fluctuations. Both materials that we have 
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discussed at length are designed ab initio, no exper- 
imental information on them yet exists to the best 
of our knowledge. 
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